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Abstract 

We reexamined the gravitational time delay of light, allowing for various models of modified 
gravity. We clarify the dependence of the time delay (and induced frequency shift) on modified 
gravity models and investigate how to distinguish those models, when light propagates in static 
spherically symmetric spacetimes. Thus experiments by radio signal from spacecrafts at very 
different distances from Sun and future space-borne laser interferometric detectors could be a 
probe of modified gravity in the solar system. 
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The nature of dark energy and dark matter has become a central issue in modern cos- 
mology. Recent observations such as the magnitude-redshift relation of type la supernovae 
(SNIa) [1] and the cosmic microwave background (CMB) anisotropy by WMAP [2] strongly 
suggest a certain modification, in whatever form, in the standard cosmological model. We are 
forced to add a new component into the energy-momentum tensor in the Einstein equation 
or modify the theory of general relativity itself [3]. Indeed, there have been a lot of pro- 
posals motivated by, for instance, scalar tensor theories, string theories, higher dimensional 
scenarios and quantum gravity (For recent reviews of modified gravity models inspired by 
the dark energy observation, e.g., [4]). Therefore, it is of great importance to observationally 
test these models. 

The theory of general relativity has passed "classical" tests, such as the dcficction of light, 
the perihelion shift of Mercury and the Shapiro time delay, and also a systematic test using 
the remarkable binary pulsar "PSR 1913-1-16" and several binary pulsars now known [5]. In 
the twentieth century, these tests proved that the Einstein's theory is correct with a similar 
accuracy of 0.1%. 

Since the time delay effect along a light path in the gravitational field was first noticed in 
1964 by Shapiro [6], this effect has successfully tested the Einstein's theory [7]. A significant 
improvement was reported in 2003 from Doppler tracking of the Cassini spacecraft on its 
way to the Saturn, with 7 — 1 = (2.1 ± 2.3) x 10~^ [8]. Here, 7 is one of parameters in the 
parameterized post- Newtonian (PPN) formulation of gravity [5]. The bending and delay 
of photons by the curvature of spacetime produced by any mass are proportional to 7 -|- 1, 
where 7 is unity in general relativity but zero in the Newtonian theory, and the quantity 
7 — 1 is thus considered as a measure of a deviation from general relativity. The sensitivity in 
the Cassini experiment approaches the level at which, theoretically, deviations 10~^ — 10~^ 
are expected in some cosmological models [9, 10]. Therefore, it is important to investigate 
the Shapiro time delay with such a high accuracy. 

In addition to the above theoretical motivation, there are advances in technologies con- 
cerning the high precision measurement of time and frequency such as optical lattice clocks 
[11] and attoseconds (10^^^ s) laser technologies [12]. ASTROD project with three space- 
crafts aims at measuring 7 at the level of 10~^ [13]. 

The purpose of this paper is to clarify the dependence of the time delay (and induced 
frequency shift) on modified gravity models and investigate how to distinguish those models 
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by using the Shapiro time delay. An important point in this paper is that we allow for 
various modified gravity theories beyond the scope of the PPN formulation. Introducing 
a new energy or length scale (e.g. extra dimension scale) may make changes in functional 
forms of the gravitational field. Thus it is worthwhile to investigate how to probe such a 
modified functional form, by using the light propagation in the solar system. Throughout 
this paper, we take the units of G = c = 1. 

In this paper, we assume that the electromagnetic fields propagate in four-dimensional 
spacetimes (even if the whole spacetime is higher dimensional). Thus photon paths follow 
null geodesies (as the geometrical optics approximation of Maxwell equation). 

Wc shall consider a static spherically symmetric spacetime, in which light propagates, 
expressed as 

ds''^-A{r)dt^ + B{r)dr^ + r^dn^, (1) 

where r and dD,'^ denote the circumference radius and the metric of the unit 2-sphere, 

respectively. The functions A{r) and B{r) depend on gravity theories. 
The time lapse along a photon path is obtained as 



, dr Bir) 1 



'0 



where h and Tq denote the impact parameter and the closest point, respectively. Their 
relation is 6^ = rQ/A(ro). 

According to a concordance between solar-system experiments and the theory of gen- 
eral relativity, we can assume that the spacetime is expressed as the Schwarzschild metric 
(rigorously speaking, its weak field approximation) with a small perturbation induced by 
modified gravity. For practical calculations, we keep only the leading term at a few AU in 
the corrections. Namely, A{r) and B{r) are approximated as 

A{r) ^ — + A^r"^, (3) 
r 

B{r) « l + ^ + S„r", (4) 

where M denotes the mass of the central body. Here, Am, Bn, m and n rely on a theory 
which we wish to test. For simplicity, we assume m — n > 0, which corresponds to a wide 
class of theories of gravity. 
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Examples of modified gravity theories are as follows. (1) n — 1/2, An — —2Bn — 
±2-\/M/r^ for DGP model with Tc is the extra scale within which gravity becomes five di- 
mensional [14]. (2) n = 3/2, An = (2/3)m^^2M/13 and Bn = -m\sj2M jl'i with graviton 
mass mg for one of massive gravity models [15, 16]. (3) n = 2, An = —Bn = —A/3 for the 
Schwarzschild-de Sitter spacetimc, that is, general relativity with the cosmological constant 
A as a possible candidate for the dark energy, though this is not a manifest modification of 
gravity. The solar system experiments are not sensitive to this model with A ~ 10~^^m~^ 
[17]. Here, it should be noted that the examples (1) and (2) give conformally fiat spacetimes 
(in the weak field approximation) and their conformal factors generate the gravitational 
time delay (and induced frequency shift), though the null geodesic in any conformally flat 
spacetime is mapped into that in the Minkowski one. 

The Cassini experiment has put the tightest constraint on the solar gravity, especially 
near the solar surface with the accuracy of 10~^ [8]. This imphes that deviations in ^(r©) 
and B{rQ) must be less than 10"^ x 2MQ/rQ ~ 10"^°, that is, \Amr^\, |S„rg| < 10"^°. 

We consider the round-trip time between pulse transmission and echo reception, denoted 
by AT. The pulse is emitted from Earth at rg, and reflected at tr. 

Up to the linear order in M, An and Bn, AT is expressed as 



AT = 2(Vr|-rg + V4 

+2M I 21n^^+v'^^ + 21n^^+^^^^ + ,/!:^+ ^^^^ 



+6t. (5) 
The extra time delay induced by a correction to general relativity is expressed as 

x(-4»^7p3T)37^ + ^"71f=Tj- W 

where we define nondimensional radial coordinates as R = r/ro. Re = r^/ro and Rr = 
t^r/t^o- For a radar tracking of a spacecraft such as Cassini, and are of the order of 
1 AU (~ 10^ km), and Tq is several times of the solar radius {Rq ~ 10^ km). Equation (6) 
can be rewritten by using special functions, though it seems less informative. Therefore, we 
take expansions of Eq. (6) in ro because of r^, ^ ro. For n 7^ 1, we obtain 

Bn An / „ I 1 „ I 1 \ Bn ~\~ A^ 
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whereas the second term of R.H.S. becomes (S„ + An) \n{rErR/rQ)rQ/2 for n — 1. 

It is convenient to use the relative change in the frequency, which is caused by the 
gravitational time delay [18], because the Doppler shift due to the receiver's motion has no 
effect owing to the cancellation at both the receipt and emission of radio signal [18]. This 
frequency shift is defined as y = —d{AT)/dt. Indeed, the frequency shift was used by the 
Cassini experiment. For a case olh <^ te-iTr, which is valid for the Cassini experiment, the 
general relativistic contribution is expressed as [5] 

Vcu = 4--. (8) 

We pay attention to the extra contribution due to modified gravity. For n 7^ 1, the extra 
frequency shift becomes 

h - -^^{rl-' + rn' -in + ^K'^yb^ (9) 

while we obtain Sy — — (A„ + S„)[ln(r£;ri?/ro) — l\bdb/dt for n — 1. Here we used dvE/dt, 
dru/dt <S dro/dt (~ db/dt) near the solar conjunction (6 <S rE,rR). The total frequency 
shift y is the sum, yoR + Sy. The impact parameter of light path changes with time, because 
of the motion of the emitter and receiver with respect to Sun. For simplicity, we assume 
that they move at constant velocity during short-time observations. The impact parameter 
changes as b{t) — ^/^'+vH^, where bo denotes the minimum of the impact parameter near 
the solar conjunction ai, t — 0, and v is the velocity component perpendicular to the line of 
sight. 

Here, we make an order-of- magnitude estimate of the frequency shift. First, we obtain 
Ugr ~ lO~^(M/M0)(r0/6)(6/f£;), where the dot denotes the time derivative, and ve is the 
orbital velocity of Earth (~ 30 km/s). The Cassini experiment reported y at the level of 
10~^^ by careful processing of the frequency fiuctuations largely due to the solar corona 
and the Earth's troposphere [8]. Multi-band measurements are preferred in order to avoid 
the astrophysical effect of the corona and interplanetary plasma on the delay, which is 
proportional to the square inverse of the frequency. 

For a receiver at tr > te, the extra frequency shift is 

,n-i7/10AUy/ (^.. + i;„)rg, \ / y-^ ( b \ fdb/dt\ 



where lOAU/ro ~ 2 x 10^. The larger the index of n, the longer the delay 5y. 

Figure 1 shows that an extra distortion due to 5y would appear especially in the tail 
parts oi y — t curves. According to the fact that no deviation from general relativity has 
been reported by the Cassini experiment [8], wc can put a constraint as 6y < 10^^^ at 
rji = 8.43AU. On the other hand, Eq. (10) gives 5y ~ 10"^^ forn = 2 and (A„ + 5„)r^ = 
O(10~^°), for instance, which are thus rejected. One can distinguish modified gravity models, 
which are characterized by various values oi n, An, Bn, from observations using receivers at 
very different distances from Sun, as shown by Fig. 1. 

Figure 2 shows the dependence of 5y on n and An + Bn- Hence, one can put a constraint 
on n and An + from dy observed. 

Equation (8) shows that the frequency shift depends only on the impact parameter b 
but not the locations of the emitter and receiver. Strictly speaking, yan still has weak 
dependence on and as shown by Eq. (5). On the other hand, 5y depends strongly on 
re and tr. The dependence of yoR and Sy on rs and tr plays a crucial role in constraining 
(or detecting) a correction to general relativity in the solar system. 

Let us imagine that time delays (or induced frequency shifts) are measured along two 
light trajectories, whose impact parameters arc denoted as 6i and 62, respectively. Then, 
we make a comparison of the two time delays. If they are in good agreement after taking 
account of a difference in the impact parameters, general relativity can be verified again. 
Otherwise, a certain modification could be required for the solar gravitational field. At this 
stage, however, one can say nothing about functional forms of the correction because the 
parameters of both n and An + Bn, which we wish to determine, enter the frequency shift. 

In order to break this degeneracy, therefore, we consider three light paths, for which 
the impact parameters of the photon paths are almost the same (several times of the solar 
radius) for convenience sake. The locations of the receivers are denoted as tri, rR2 and trs, 
where the subscripts from 1 to 3 denote each light path. We assume that is constant in 
time for simplicity. It is a straightforward task to take account of the eccentricity of the 
Earth orbit and a difference between the impact parameters. 

We make use of a difference such as 1/2 — yi and 7/3 — yi, in order to cancel out general 
relativistic parts. We find 



Z/2 -yi 



An + B, 
n-1 



n 



ir 



.n-l 
Rl 



rR2')b 



dh 
di' 



(11) 
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It should be noted that 1/2 — yi is proportional to An+Bn- Hence, the following ratio depends 
only on n as 

1/2-1/1 r-^-r^^^'' ^ ' 

Thereby, one can determine the index n. Next, one obtains A„ + Bn by substituting the 

determined n into Eq. (11). 

In summary, we have clarified the dependence of the gravitational time delay on modified 
gravity models. For neighboring light rays, Eq. (7) gives almost the same value so that one 
can hardly distinguish models of gravity. This implies that we should prepare receivers at 
very different distances from Sun. Our result could be used for exterior planets explorers 
such as New Horizons, which were launched in 2006 and their primary target is Pluto and 
its moon, Charon at distance from Sun ~ 40 AU [19]. In future practical data analyses, 
however, it would be safer to use the original integral form as Eq. (6), because Eq. (7) is 
an approximate expression. 

Furthermore, b becomes the same order of r^, r^; for future space-borne laser intcrfcromct- 
ric detectors such as LISA [20], DECIGO [21] and especially ASTROD [13]. These detectors 
are in motion in our solar system. Namely, tr and b change with time. Therefore, the so- 
phisticated experiments by space-borne laser interferometric detectors, which are originally 
designed to detect time-dependent part of gravity, i.e. gravitational waves, could probe also 
a time-independent part of gravity at the relative level of Ay ~ Au/u ~ AL/L < 10~^°. It 
would be important to make a feasibility estimate for these detectors. Clearly, a stronger 
test can be done by not a single experiment but combining several ones. In addition, more 
precise measurements of the Shapiro time delay with binary pulsars may put a constraint 
on the modifications discussed in this paper, especially in the strong self-gravitating regime. 
Further investigations along these lines will be done in the future. 
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FIG. 1: Dependence of the frequency shift on the distance vr and the index n. The long dashed, 
short dashed and dotted curves denote the frequency shift for {n,rR) = (3/2,10AU), (n,rR) = 
(2, lOAU), {n,rji) = (2, lAU), respectively. The long dashed curve for n = 3/2 and rR = 10 
AU is overlapped with the solid curve denoting the general relativistic case. Here, we assume 
(A„ + B„)rg = 3x 10-11. 
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FIG. 2: Contours of 5y on the n - |A„ + Bn\rQ plane. The solid, long-dashed and short-dashed 
curves correspond to Sy = 10~^^, 10~^^, 10~^°, respectively, where we assume = 1 AU, vr = 40 
AU, h Tq and dh/dt ^ ve- The limit due to the current technology is 5y ~ lO"-*^^. The shaded 
region above the dotted curve {5y = 10~^^ for tr = 8.43 AU) is excluded, because no devitation 

up to 0(10^^^) has been detected by the Cassini experiment [8]. 
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